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ABSTRACT 

The evolution of forces acting on horizontal cylinders subjected to 
impact by an harmonically oscillating free surface has been investigated 
both theoretically and experimentally. 

The experiments were carried out in a large U-shaped water tunnel, 
with cylinder diameters ranging from 3.0 to 6.5 inches. The results 
have been expressed in terms of two force coefficients. The first is 
the slamming coefficient which expresses the normalized force acting on 
the cylinder at the time of impact; and the second is the maximum drag 
coefficient which occurs when the cylinder is immersed approximately 
1.75 diameters in water. The slamming coefficient was found to be equal 
to tt, as predicted theoretically. It was also found that the slamming 
coefficient may be amplified to a value as high as 2rr through the 
dynamic response of the cylinder. The maximum of the normalized force 
in the drag dominated region was found to be approximately equal to 2.0. 
It decreased with increasing Froude numbers to a value of about unity. 

It is recommended that further experiments be conducted with sand- 
roughened cylinders and the dependence of the rate of rise of the impact 
force on the flow parameters be investigated. 
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NOMENCLATURE 



A Amplitude of the wave motion 

A. Immersed cylinder area 

A. Normalized immersed cylinder area 

Buoyancy corrected drag coefficient 
C s Slamming coefficient 

C° Initial value of the slamming coefficient 

D Cylinder diameter 

F Instantaneous total force 

Fr Froude number 

f Natural frequency (Hz), oo = 2ir f 

g Gravitational acceleration 

K Keu I egan-Carpenter number 

L Length of cylinder 

M Mass 

m Added mass 

m Normalized added mass 

p Pressure 

Re Reynolds number 

r Rad i us of cy I i nder 

s Distance between pressure taps 

T Wave period 

t Time 

t Rise time 
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Maximum free surface velocity 

Instantaneous velocity 

Initial velocity 

Depth of cylinder immersion 

Normalized depth of cylinder immersion 

Damping coefficient 

Instantaneous height of wave surface above 
mean water level 

Height of cylinder with respect to mean water level 

Fluid density 

Velocity potential 

Natural frequency (rad/sec) 

V i scos ity 
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I . I NTRODUCT I ON 



Information concerning the forces acting on bluff bodies subjected 
to wave slamming is of significant importance in Mechanical Engineering, 
Ocean Engineering, and Naval Arch i tecture. The design of structures 
which must survive in a wave environment is dependent on a knowledge of 
the forces which occur at impact, as wel I as on the dynamic response of 
the system. Two typical examples include the structural members of off- 
shore drilling platforms at the splash zone and the often-encountered 
slamming motion of ships. 

The general problem of hydrodynamic impact has been studied exten- 
sively Cl 3 motivated in part, no doubt, by its importance in ordnance 
and missile technology. Extensive mathematical models have been developed 
for cases of simple geometry such as spheres and wedges, and these models 
have been well supported by experiment. Unfortunately, the special case 
of wave impact has not been studied as extensively. Kaplan and Silbert 
[ 22 developed a solution for the forces acting on a cylinder from the in- 
stant of impact to full immersion. Dalton and Nash [33 conducted slam- 
ming experiments with a 0.5 inch diameter cy I inder with sma I I amp I iTude 
waves generated in a laboratory tank. Their data exhibited large scatter 
and showed no particular correlation with either the predictions of the 
hydrodynamic theory or identifiable wave parameters. It appears that the 
evaluation of the slamming effects with wavy flows is extremely difficult 
partly because of the limited range of wave amplitudes that can be achieved 
and partly because of the difficulty of measuring the fluid velocities 
at the instant of impact. 
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In essence, the present research was undertaken to extend the work of 
Dalton and Nash and to improve the experimental method by using a large 
U-shaped water tunnel. In addition, it was intended to demonstrate a cor- 
relation between theory and experiment by taking into account the dynamic 
response of the system. Specifically, the goals of the project were as 
fol lows: 

A. to examine the existing theoretical models in determining wave 
slam forces on circular cylinders; 

B. to furnish data, obtained under controlled laboratory conditions, 
about the forces acting on smooth circular cylinders subjected to 
impact with an harmonically oscillating water surface; 

C. to determine the relative importance of the inertia dominated and 
drag dominated forces during wave impact; and 

D. to correlate such data with respect to identifiable wave parameters 
such as the Froude number, Fr; the Keu I egan-Carpenter number, K; 
and the Reynolds number. Re. 
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II. THEORETICAL ANALYSIS 



The traditional approach used in the design of off-shore structures 
involves the classical Morison equation to determine the forces due to 
wave motion. However, wave impact is generally defined as the early 
stages of the penetration of a sol id body into a wave surface. At the 
instant of contact the fluid in the vicinity of the body undergoes large 
accelerations which give rise to large forces. As the body becomes more 
fully immersed, forces due to viscous drag and separation effects become 
predominant. Thus the inertia and drag coefficients used in the Morison 
equation are not constant, and the problem becomes very difficult, even 
for simple geometries. 

The general case of hydrodynamic impact is usually described by using 
incompressible potential flow theory. For the case of a moving body with 
mass, M, and velocity, v q , impacting a quiescent surface, the system momentum 
is Mv q . Neglecting non-conservative forces, the momentum of the system 
is unchanged during penetration. However, the mass of the system 
increases due to the fluid which is set in motion in the vicinity of the 
body. Also known as "added mass", m, it results in reduction of the 
velocity. Thus the system momentum after penetration is (M + m) v = Mv q . 
The impact force at any instant is a function of m and 3m/3t. The 
so I ut ion, therefore, requires knowledge of the added mass and its time 
der i vat i ve. 

It should be noted that viscous effects, bei ng non-conservative, will 
alter this somewhat. Additionally, high speed entry (as with projectiles) 
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may require that compressibility effects of both the fluid and the air 

above it be taken into account. However, for low speed entry, these 

effects are generally neg I ected, and the fluid motion is described by a 

2 

velocity potential^ , which satisfies V <j) .= 0. The added mass term 
can then be calculated. The problem remains extremely difficult, though, 
because the free surface tends to "pile up" around the body. This 
results in an unsteady flow problem with a time dependent boundary condi- 
tion (the free surface). Approximate solutions have been obtained £4] 
for simple shapes (spheres, cones, cylinders, wedges, flat plates, etc.) 
by assuming a constant plane free surface and then applying a correction 
for the distortion. 

Kaplan and S i I bert C2H developed the following relation for the force 
acting on a horizontal circular cylinder by a wave system which propa- 
gates normal to it, (See Fig. |.); 

■£■ = pg a. + (m + p Aj ) n + n 2 ( I ) 

in which F represents the force acting on the cylinder; L, the length 
of the cylinder; p, the fluid density; g, the gravitational acceleration; 
A., the immersed area; m, the added mass per unit cylinder length; q, the 
instantaneous height of the wave surface above the mean water level; and 
z, the instantaneous depth of cylinder immersion. The first and second 
derivatives of q with respect to the time are denoted by q and q. The 
added mass is given by Taylor £5] as: 

m = j pr L-j-( 2 YV - Q - > + I -cos 8) + (sin8-9)J (2) 

in which r represents the radius of the cylinder, and 8 is defined as 
shown in Fig. I. 
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Fig. I. Geometry for theoretical analysis of wave impact. 
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The motion of the free surface, q, is related to the maximum ampli- 



tude by: 



r| - A sin 2irt/T 



(3) 



where A and T represent the amplitude and period of the free surface 

oscillations. Equation (I) can also be written in the form of a slamming 

2 

coefficient, C = 2F/pDL U as: 

' s ^ m 



„ t gr , r . r . 2tt , , 3m 2 2tt , 

C = A . =• -(m + A.)-7- s| n-^r-++ — cos t 

s i ,, 2 i A T T 



(4) 



U 



m 



3z 



where 



2 2 
A. = A./r , m = m/p r 
ii 

z = z/r , U = 2 ttA/T 
' m 

With z = r(l-cos -^) , Equation (2) can be used to relate the added 
mass to the depth of immersion. Figure 2 is a plot of 3m/3z as a func- 
tion of z/D. Clearly, 3m/3z begins with an initial value of tt and drops 

2 

rapidly. The quantity gr/U m is related to the Froude number, Fr, by 

Fr = U 2 /2gr. Thus C = f(Fr, £ , — ) . 

3z 

The rate of change of the normalized added mass with z depends on 0 
and hence on the time measured from the instant of impact (See Fig. 2). 
For example, for the case where ri Q = 0 and t = 0, 3m/3z = tt and C° = tt. 
Consequently, the for the particular case under consideration, C° at the 
instant of impact does not depend on either the size of the cylinder or 
the flow parameters. 
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Fig. 2. 3m/3z as a function of z/D. 
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The fact that the cylindrical members of a structure in the splash 
zone are not necessarily located just at the mean water level requires 
the determination of the particular value of r| for which the slamming 
force at the time of impact is a maximum. It can be demonstrated through 
the use of Equation (4) that the maximum impact force occurs for the 
case of iq = 0. For this purpose Equation (4) was evaluated with the aid 
of a computer program which allowed variations in r/A and h Q /A. C g was 
then plotted as a function of z/D from zero to unity. Figure 3 is an 
example of one of these plots and is typical for most cases. The variables 
were: 



r/A = 0. 125 

n Q /A = 0.0, 0.4, and 0.8 

As can be seen readily from Fig. 3, C s is largest at z/D = 0 for lq^/A =0.0 
and starts at a value of tt and drops to a minimum at z/D of approximately 
0.5. It rises again as z/D approaches 1.0. Thus Equation (4) indi- 

cates that C , and consequently the impact force, is of an impulsive 
nature beginning with a finite value (not zero) at the instant of impact. 
Since viscous effects are neglected, one would expect the solution to 
deviate from the actual situation as the cylinder becomes more fully 
immersed. Where this becomes the case can only be determined by experiment. 

However, one can assume that the solution would be valid for small 
depths of immersion, which is the case for the instant of impact up to 

z/D of 0.02 or so. With this restriction, A. and m are sma I I, as is 

2 

sin 2irt/T. Additionally, cos 2iTt/T is very nearly equal to 1.0. There- 
fore C reduces to: 
s 

C ~ C° = (3m/3z), _ 
s s t=0 
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Fig. 3. IheoreTical values of C_ as a -unction of 

z/D for several values of n /A. 

o 
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It is not realistic to assume that the impact force rises from zero 
to ir instantaneously. Several factors, specifically the compressibility 
of the air between the cylinder and water surface, entrapped gases in 
the water, and surface i rregu lari ties would account for some finite rise 
time. Nonetheless, the rise time can be expected to be short, i.e. in 
the order of milliseconds. The exact nature of the rise is an interest- 
ing question for further study. However, in this analysis C s is assumed 
to vary I inearly during a rise period, t r . For t _> t , C s drops from ir 
as 3m/3z. Figure 4 is a representation of this assumption. Exactly how 
lonq of a time interval t is will be discussed later. 

The realization that the impact force is of an impulsive nature requires 
consideration of the fact that this force does not act on a perfectly rigid 
body, but rather on a cylinder which is supported elastically. The re- 
sponse of such a system approaches that of a rigid body only if its 
natural frequency approaches infinity. Additionally, the response of the 
system to an impulsive force is heavily dependent on the exact nature of 
the force itself as well as on the system natural frequency. 

In general, the instantaneous d i sp I acement, x(t ), of a system of mass, 
M,with a spring constant, k, i s given by (See e.g. C6]) : 



in which £ represents a dummy variable; F(t), the driving force; and 
g(t), the response to a unit step excitation. Equation (6) is readily 
recognized as the DuHamel superposition integral, which can be expanded 



t 




( 6 ) 



o 



as fol lows for g(t) = — sin w t: 

a Mw n 

n 
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Fig. 4. C (+) as a function of t. 
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X ( t ) 



I 

Moo 

n 



cos O) t 
n 




o 



s i n co tdt 
n 



+ 



Mco 



s i n 



CO t 

n 




t 

F( + ) cos co tdt 
n 



( 7 ) 



o 

in which co is the natural frequency of the cylinder and supports. 

Letting F(t) = F f(t) and changing the variable of integration to 

9 = co t, one obta i ns 
n 



k x(t) 
F 

o 



9/u 



- cos 9 



/ 



f (t) si n 9d9 + sin 9 



9/co 

r " 

/ f(t)cos 9d9 
o 



( 8 ) 



in which k = co n M, the spring constant. The term on the left-hand side 
can be interpreted as the ratio of the instantaneous force sensed by the 
supports of the cylinder to the actual mean force acting on the cylinder. 
In other words, the response of the system at time t is k x(t) which may 
be different from F q . For the simple case in which F(t) is a step 
function as shown below. 



F 

o 



t 



equation (8) reduces to: 



k x(t) 
F 

o 



( I - COS CO t ) 

n 



( 9 ) 
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It is apparent from Equation (9) that the instantaneous force sensed by 
the system can be anywhere from zero to two times the actual mean force. 

If one assumes, as before, that the impact force is as described by 
Fig. 4, then Equation (8) must be evaluated using C s (t) as F(t). Addi- 
tionally, damping can be taken into account by writing g(t) as C63: 

g(t) = e ^n" 1 " sin w t (10) 

n 

where £ is the damping coefficient. Thus Equation (8) can be rewritten 
by replacing the forcing function, F(t), by C s (0) as: 



0/w 



k x(t) 
C ° 



= - COS0 



/ 



0/(0 



C (0)e ^ 0 sin0d0 t cos0 



/ 



C s (0)e"^ 9 cos0d0 



(II) 



Equation (II) was solved with the aid of a simple computer program which 

used a trapezoidal integration of variable step. Values of (o n were 

taken as 358 sec ' and 628 sec ' which corresponded to the measured values 

of (o for a 6.0 inch and a 3.0 inch diameter cylinder respectively, as 

discussed later. The damping coefficient was 0.014 in both cases, also 

correspond i ng to measured values. The rise time was varied from zero to 

approximately 0.025 seconds. Figures 5 through 8 are representative plots 
k x ( i~ ) 

of for the 6 inch cylinder, for various values of t . 

C ° 

u s 

Figure 5 represents Equation (II) plotted for a rise time of 0.0001 

second; Fig. 6, for a rise time of 0.0100 second; Fig. 7, for a rise 

time of 0.0195 second; and Fig. 8, for a rise time of 0.0230 second. 

k x ( "t" ) 

Note that for a very short rise time the value of reaches a value 

C ° 

5 

of approximately 1.7 at t = 0.001 seconds and then drops off. As the rise 
time increases, the value of the first peak drops off, with a "double 
peak" appearing at a rise time of 0.0195 seconds. 
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for t = 0.0001 second 
r 
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Fig. 6. kx(t)/C° as a function of t 

for t = 0.010 second, 
r 
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Fig. 7. kx(t)/C° as a function of t 

s 

f° r c r = 0.0195 second. 
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Fig . 8 . kx ( t Vc° as a function of t 

s 

for t = 0.023 second, 
r 



26 



The interpretation of the above results indicates that, depending on 
the rise time of C s (t), values of the apparent slamming coefficient, C°, 
may lie between roughly 0.5 and 1.5 of the theoretical instantaneous 
value of tt. Again this applies only to a small depth of immersion, cor- 
responding to the initial moments of impact. The significance of this 
fact is that values of C° measured in the laboratory may show wide scat- 
ter depending on the rise time for any given experiment. Additionally, 
if the surface is not perfectly plane, rise times may vary from experiment 
to experiment resulting in an apparent non-repeatability. 

Figure 9 is a plot of k x(t)/C° as a function of f t for both the 
first and the second peak, using co n values for both a 6.0 inch and a 3.0 
inch cylinder. From this figure, the drop in k x(t)/C° for the first 
peak with increasing rise time is apparent. It is also evident that the 
effect of to is confined to a narrow range of f t values smaller than 
about 0.5. The force acting on the cylinder is amplified by the dynamic 
response of the system for a value of f t smaller than 0.9 and 
attenuated for f t^ values larger than 0.9. The second peak occurs 
only for f t values larger than 1.09. The value of k x(t)/C° increases 

from 0.84 to about unity as f t increases from 1.09 to 1.5. 

7 n r 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 



The equipment used to generate the harmonically oscillating flow has 
been used extensively at this facility over the past two years. The 
apparatus is described in reference C7j. The salient features, 
as well as the adaptation for this work, are briefly described in the 
fo I lowing. 

A. U TUNNEL 

Figure 10 is a schematic diagram of the U tunnel. A butterfly valve 
arrangement at the top of the left leg allows that side to be completely 
sealed off. Compressed air can then be admitted until the water level 
is driven to the desired height in the opposite leg. Rapid opening of 
the butterfly valve allows the water to oscillate at a natural frequency 
determined by the geometry of the tunnel. The tunnel contains approxi- 
mately 5000 gallons of water and the period, T, is 5.5 seconds. Ampli- 
tudes up to 4.0 ft can be generated. 

Previous work with this tunnel has demonstrated its capability in 
studies involving oscillating flow. Measurements have confirmed that 
separation does not occur at the corners and there are no discontinuities 
where the individual sections are joined, which might otherwise lead to 
disturbances in the flow. These factors were important to this project 
since uniformity of the free surface was a major concern. 

B. TEST CYLINDERS 

Three, five, six, and six and one half inch diameter aluminum 
cylinders were used in measuring impact forces. Measuring 2.98 ft in 
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Butterfly Valve 
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Schematic diagram of the U tunnel. 




length, they were constructed of aluminum pipe, turned to the final 
diameter, and polished to assure a hydrodynamica I I y smooth surface. 

Caps fitted to each end prevented water leakage into the cylinder. 

Each cap contained a double precision ball-bearing mounted flush with 
the face. As will be noted later, the force transducers were attached 
to the cylinder via these bearings, and the purpose was to allow free 
rotation of the cylinder. 

C. FORCE MEASUREMENTS 

Two cantilever-beam force transducers were used to measure the in- 
stantaneous in-line (drag) and transverse (lift) forces. The output 
consisted of an electrical signal from strain gages attached to the 
beam. The rounded end of the transducer assembly was designed to fit 
snugly in the bearings in the ends of the cylinders. Neoprene "0" 
rings were added to the end of the transducers to provide axial position- 
ing of the cylinder within the tunnel. 

Calibration of the transducers was accomplished by hanging a load 
from the center of the cylinder. This not only established the level 
of the electrical output for a known load but also reaffirmed that the 
transducers were correctly oriented. 

As with other elements of the experimental apparatus, these trans- 
ducers had been used frequently at this facility for over two years and 
the description of their construction may be found in reference [73. 

No changes have been noted in the calibration of these transducers since 
their installation, a fact which establishes some degree of confidence' 
in their use. 

Initial experimental efforts established that the transverse forces 
were very small compared to the in-line forces, and consequently measure- 
ment of these forces was discontinued early in the study. 
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D. MEASUREMENT OF FLUID VELOCITY 



Two pressure taps, one in each leg of the tunnel, were connected by 
tubing to a differential pressure transducer located midway between the 
taps. The acceleration of the fluid can be determined from Ap = ps ^pp , 
in which Ap represents the pressure differential; p, the fluid density; 
s, the distance between the taps; and — , the instantaneous acceleration 
of the fluid. The maximum velocity and amplitude are related to the 
acceleration by: 



dt m 



2rr 

T 



U = & 2 

m T 



( 12 ) 



E. RECORDING EQUIPMENT 

A three channel strip chart recorder was used to record simultaneously 
the output of the in-line force transducer and the differential pressure 
transducer. The third channel of the recorder was also used to record 
the force measurements, but the signal was first passed through a low- 
pass filter to remove the oscillation of the cylinder at its natural fre- 
quency, thus providing an "average" force record. 

Inasmuch as the impact force to be measured consisted of a rather 
short rise time, the bandwidth of the recorder was of some importance. 

The frequency response of this recorder was f I at to 100 Hz. In order to 
verify the response under the test conditions, several experimental runs 
were made in which the force was both recorded using the strip chart re- 
corder and a tape recorder simultaneously. The taped signal was then 
fed back to the strip chart recorder at 1/4 speed. This method assured 
that the frequency response of the strip chart recorder was adequate. 



32 



F. EXPERIMENTAL PROCEDURE 



The success or failure of any experimental work depends on careful 
attention to those factors which may have a significant effect on the 
measurements being taken. For this reason, extreme care was taken with 
all flow parameters. 

Prior to filling the tunnel, the cylinder was struck lightly, and the 

natural frequency and damping coefficient were determined. The tunnel 

was then filled, immersing the cylinder several diameters, and again the 

natural frequency was recorded. Thus the value of co n was obtained both 

in air and in water. In the calculations the value of co in air was 

n 

used since at impact and shortly thereafter, the cylinder behaves as 
though it were not immersed. 

Since the height of the mean water level with respect to the bottom 
of the cylinder was shown by the analysis to be important, the water 
level in the tunnel was adjusted carefully by slowly filling or emptying 
the tunnel until a slight ripple occurred at the water surface due to 
contact with'the cylinder. This method assured that the mean water level 
coincided with the cylinder bottom surface. 

Of paramount importance was the condition of the free surface at 
the instant of impact with the cylinder. Ideally the surface should be 
perfectly plane with no surface disturbances. It is this requirement 
that lends a certain degree of difficulty in the experimental procedure. 
Until the fluid in the tunnel has completed one cycle or so, the period 
of oscillation is not established. This requ i res that the cylinder be 
immersed at least once before the period of oscillation has settled down, 
resulting in some surface disturbance. The observations of the wave 
impact indicated that one slam generated minimal surface disturbances, 
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with the disturbances increasing with subsequent impacts. For this 
reason, all data were recorded for the first impact only, and successive 
impacts were discounted. Approximately ten minutes were allowed between 
each experimental run in order to a I low the water surface to reach a 
quiescent state. 

The stripchart recorder was run at a speed of 200 mm per second, each 
division representi ng 0.005 second. Figure II is a typical example of 
the data obtained. The top trace is the filtered force trace, and the 
bottom one is the unfiltered force trace. Figures 12 and 13 show addi- 
tional unfiltered force traces for various initial flow conditions. 



G. REDUCTION OF DATA 

Two force coefficients were defined in this investigation. They are: 



C 

s 



2F 

pDLU 2 
m 



(13) 



and 

r = — - 

h pDLU 2 2U 2 
r m m 



(14) 



The first of these refers to the slamming coefficient as defined by 
Equation (4). The second coefficient represents the second maximum of 
the normalized force with the buoyancy subtracted. This maximum occurs 
after the cylinder is fully immersed. 

As mentioned previously, the stripchart recorder was run at 200 mm 
per second, thus each mm corresponded to 0.005 second. The measured 
force was then read each 0.005 second up to 0.5 seconds or so, depending 
on the diameter of the cylinder. These points were then evaluated with 



34 




35 



Fig. 11. Data trace showing filtered and unfiltered traces. 
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Fig. 12. Data trace for a run with double first peak. 
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Fig. 13. Data trace for a run with decreased first peak. 



the aid of a computer program which solved Equation (4) and plotted 
C s as a function of z/D. Figure 14 is typical of such plots, which shows 
the measured value of C , the theoretical value of C s from Equation (4), 
and the normalized buoyant force. 

The filtered force trace obscured the magnitude of the initial impact. 
However, after the initial peak, the filtered trace was found to corre- 
spond to the average value of the unfiltered trace. Therefore, the data 
points during the early stage of impact were read from the unfiltered 
trace, and the remainder were read from the filtered trace. The main rea- 
son for this is the obvious difficulty in interpreting the average value 
which is obscured by the oscillation of the cylinder. 
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Fig. 14. C s (theoretical), (experimental), and buoyancy as 
function of z/D. 



IV. DISCUSSION OF RESULTS 



The initial value of the slamming coefficient, C 



,o 




s 



function of various parameters such as the Froude number, Fr = 



Reynolds number. Re = U m D/v ; and the Keu I egan-Carpenter number, K = 

U T/D. These plots have revealed that there is no identifiable corre- 
m 



tion of this lack of correlation that led to the consideration of the rise 



i ng f t values were determined from Fig. 9. Figure 15 is a plot of 



there is considerable scatter, this figure shows that for small values 
of U , or large values of D, the rise time is larger and thus yields 
smaller amplification (see Fig. 9.). For larger impact velocities or 
for smaller cylinders, the rise time is considerably smaller, and the 
amplification factor approaches 1.7. If should be noted 
once more that surface i rregu I ar i t i es , dissolved air, and other experi- 
mental uncerta inties do not permit a more accurate determination of the 
rise time. One must add that the rise time may vary also with the 
physical character i sti cs of the fluid above the free surface, the angle 
of inclination of the cylinder, the length to diameter ratio, etc. Be 
that as it may, the Froude number appears to be the dominant parameter 
as in all free surface flow problems. 

The physical reasoning behind the correlation of the normalized rise 
2 

time, f t , with gD/U is as follows. The disturbances in the free 
nr 3 m 



I at ion between C° and the said parameters. In fact, it is the rea I i za- 



time, t^. Subsequently, the amplification or attenuation of C° was calcu- 
lated by dividing C° by its theoretical value of it. Then the correspond- 



2 

f t as a function of the inverse Froude number gD/U . Even though 
nr 3 m 3 
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surface during the first passage of the flow increase with increasing 

diameter. On the other hand, the larger the amplitude, the longer the 

time for the disturbances to attenuate. Thus, the rise time decreases 

with increasing amplitude. Finally, the larger the acceleration of the 

2 

flow, i.e. A(2 tt/T) relative to the gravitational acceleration, the 

faster is the rate of attenuation of the disturbances C8,93* A simple 

dimensional analysis shows, therefore, that f t increases with 
’ nr 

2 

gD/A(T /A). Since U is proportional to A/T, one finds that f t is 

m nr 

2 

a function of gD/U 

A careful comparison of the measured and calculated force traces 
during the early stages of impact lends ample credence to the idea of 
rise time. Figures 5 through 8 and II through 13 show that for very 
short rise times there is a single sharp peak in the force trace. For 
intermediate rise times (See Figures 7 and 12.) a double peak occurs. 
Finally, for larger rise times (See Figures 8 and 13.) the first peak 
practically disappears, and the response to the slamming force is 
attenuated by the system. 

It is event from the foregoing that the determination of the initial 
slamming force is extremely difficult because of its magnification or 
attenuation by the dynamic response of the system. The evidence presented 
above shows that the normalized slamming force should be assumed equal to 
it and that the response of the system be determined through a straight- 
forward vibration analysis. Such an analysis will yield the largest 

force amplification for a normalized rise time of f t — 0.25. For 
v nr 

design purposes, the maximum slamming coefficient may be as large as I.7 tt. 

Following the initial impact, the net force acting on the cylinder 
begins to decrease since the Om/Oz contribution decreases. During 
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this stage of the flow the cylinder undergoes damped oscillations at its 
natural frequency. As the free surface rises further, the buoyant force 
i ncreases, and the separation effects begin to give rise to larger drag 
forces. The buoyancy subtracted fluid force reaches its maximum at z/D 
values from about 1.5 to 2.0. Even though the present flow situation, 
in which there is a free surface, cannot be directly compared with an 
impulsive flow about a circular cylinder, the rise of the drag force to 
a maximum at z/D— 1.75 is very much like the rise of the drag coeffi- 
cient to a maximum at a relative fluid displacement of 2.0 in impulsive 
flow [|o]. The said rise in the drag coefficient is because of the for- 
mation of two symmetrical vortices behind the cylinder. As the motion 
continues, the vortices become asymmetrical and shed a I ternat i ng I y. 

The maximum force correspond i ng to the second peak has been 

I 2 

corrected for buoyancy, normalized by tc pU^ DL,and plotted as a 
function of the Froude number in Figure 16. Two facts are apparent: 
First, the scatter in the data is considerably less than that correspond- 
ing to the initial impact; second, C h [See Equation (14).] decreases 
with increasing Froude number to a value of about unity and remains 
nearly constant for all Froude numbers from 0.5 to 2.0. 
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Fig. 16. C, as a function of Fr. 



V. CONCLUSIONS 



The theoretical and experimental investigation of flow impact on 
circular cylinders warranted the following conclusions; 

1. The dynamic response of the system is as important as the impact 
force and one cannot be determined without taking the other into 
consi deration. 

2. The initial value of the slamming coefficient is essentially equal 
to its theoretical value of tt, and the system response may be 
amplified or attenuated depending upon its dynamic characteristic. 

3. The normalized rise time appears to be a function of the Froude 
number. 

4. Following the impact, the cylinder undergoes damped oscillations 
at its natural frequency. 

5. The buoyancy corrected normalized force in the drag dominated 
region reaches a maximum at a relative fluid displacement of 
about z/D — 1.75. Subsequently, the shedding of vortices as 
well as the deceleration of flow reduces the maximum drag coeffi- 
cient to below unity. 

6. Since it does not appear that the rise time can be determined 
theoretically, the dynamic response of the system should be 
analyzed using the theoretical value of the slamming coefficient 
and the experimentally determined rise times. 
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